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To  account  lor  range  errors  obtained  when  a code-modulated  radio  signal  is  used  to  make 
range  mcasureinei  ts  in  a scattering  environment,  a mechanism  termed  “carrier  cancellation"  is 
proposed.  This  allows  for  destructive  interference  between  several  nearly  direct  path  signals, 
leaving  the  possibility  for  a weaker  signal  that  has  traveled  over  a longer  path  to  predominate 
A mathematical  formulation  of  the  phenomenon  is  given,  and  the  equations  derived  are  used 
for  a computer  simulation  of  the  effects.  The  results  of  the  simulation  confirm  that  the  pro- 
posed mechanism  represents  a viable  explanation  ol  the  observed  range  errors.  t - — 
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PREFACE 


The  work  reported  here  was  undertaken  in  connection  with  a position  deter- 
mining system  that  uses  the  ranges  from  a receiver  on  the  ground  to  several  satellites 
having  known  positions  to  determine  coordinates  of  the  ground  location.  These  ranges 
are  proportional  to  the  delay  time  between  transmission  by  the  satellite  and  reception 
on  the  ground.  The  delays  are  determined  by  correlation  with  a signal  generated  on 
the  ground  that  is  a replica  of  the  satellite  signal.  Since  it  is  desirable  to  use  the  system 
in  the  presence  of  vegetation  and  other  potential  scatterers,  it  is  important  to  know  the 
characteristics  and  extent  of  errors  that  can  be  introduced  by  multipath  signals  caused 
by  such  scatterers.  Field  measurements  had  been  performed  in  the  presence  of  vegeta- 
tion and  the  results  indicated  that  anomalous  range  measurements  sometimes  occur. 
The  present  study  is  intended  to  provide  a plausible  explanation  for  anomalous  range 
measurements  obtained  in  such  a scattering  environment. 

The  work  was  performed  in  the  Center  for  Theoretical  and  Applied  Physical 
Sciences  in  the  Research  Institute  of  the  U.S.  Army  Engineer  Topographic  Labora- 
tories, under  Project  4AI6J  J02852C,  Task  S3,  Work  Unit  175 2CS30003. 

The  authors  are  especially  grateful  to  Mr.  Jack  Jacobson  and  his  co-workers  for 
making  available  the  results  of  the  field  experiments  in  a reduced  and  usable  form. 
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SCATTERING  OK  A CODE  MODULATED  RADIO  SIGNAL 


AND  ASSOCIATED  MULTIPATH  RANGE  ERRORS 


INTRODUCTION 

Objective.  In  making  range  measurements  using  a code-modulated  radio  signal  in 
the  presence  of  scattcrers,  sometimes  large  range  measurement  errors  were  obtained 
In  all  cases,  the  range  measurements  obtained  were  larger  than  the  actual  range.  Often, 
a proper  measurement  ol  the  range  could  be  restored  if  the  receiving  antenna  was 
moved  less  than  three-tenths  of  a meter  (about  I foot).  Figures  I through  3 show 
examples  of  actual  measurements.  The  purpose  of  this  study  was  to  develop  a mathe- 
matical model  that  would  account  for  the  anomalous  range  measurements  and  would 
also  provide  for  computations  leading  to  a numerical  simulation  of  these  measure- 
ments. The  computations  were  to  be  performed  for  some  representative  cases  to  dem- 
onstrate the  plausibility  of  the  model  and  to  give  some  insight  into  the  mechanisms 
and  relative  significance  of  the  various  parameters  involved. 

The  basic  scattering  phenomena  upon  which  the  characteristics  of  the  received 
signal  ultimately  depend  were  not  investigated  because  such  an  undertaking  would 
have  exceeded  the  scope  of  this  study.  It  is  sufficient  to  know  that  diffraction  and 
scattering  of  radio  waves  in  foliage  will  generate  interference  patterns  with  spatial 
distributions  of  electromagnetic  energy  that  show  considerable  local  variations.  For 
that  reason,  the  effects  of  scattering  were  included  only  in  the  form  of  delayed  multi- 
path  signals  with  somewhat  arbitrarily  assigned  delay  times  and  amplitudes.  Also,  no 
attempt  was  made  to  model  the  receiver  in  detail,  but  it  was  considered  to  possess 
idealized  characteristics  of  receivers  in  general. 

The  Experimental  Configuration  to  be  Simulated.  The  distance  between  two 
points  is  measured  by  determining  the  travel  time  for  a signal  between  a transmitter 
and  receiver  as  shown  in  figure  4.  A biphase  code-modulated  signal  sent  by  the  trans- 
mitter is  received  and  correlated  with  an  identical  signal  produced  at  the  receiver.  By 
using  an  accurately  synchronized  system  of  clocks,  the  time  delay  is  obtained  and 
hence,  the  corresponding  range  or  distance  from  the  transmitter  to  the  receiver. 

When  such  a system  is  used  in  an  environment  containing  scatterers  such  as  trees 
and  foliage,  several  signals  transmitted  over  different  paths  can  possibly  be  received 
(figure  5).  In  such  surroundings,  it  is  important  that  the  wave  propagated  over  the 
most  direct  path  transports  sufficient  energy  that  the  receiver  may  lock  onto  this 
component;  otherwise,  the  receiver  may  lock  onto  another  component  and  give  a 
faulty  range.  One  mechanism  that  can  cause  this  to  happen  is  destructive  interference 
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Figure  2.  Second  Example  of  Range  Errors  Caused  by  Carrier  Cancellation. 


Figure  3 Third  Example  of  Rang*  Error*  Canted  by  Carrier  Cancellation. 


TRANSMITTER 


RECEIVER 


Figure  4 Range  Measurement  (or  a Single  Direct  Path. 


Figure  5.  Range  Measurement  in  the  Presence  of  Scatterers. 

between  several  of  the  more  direct  waves,  leaving  the  strongest  signal  present  to  be 
one  that  has  been  propagated  over  a longer  path  length.  Because  this  interference  is 
associated  with  the  wavelength  of  the  carrier,  it  has  been  termed  “carrier  cancellation.” 
To  explain  how  and  when  it  occurs,  an  understanding  of  the  detailed  structure  and 
properties  of  the  modulated  signal  is  necessary. 

Signal  and  Code  Characteristics.  The  transmitted  signal  consists  of  a carrier  with 
an  impressed  modulated  signal  with  a repeated  code  or  string  of  pulses  having  a value 
+ 1 or  -I  (figure  6). 


+1 


-1 


Figure  6.  Modulated  Signal. 


At  each  discontinuity  of  the  modulated  signal,  the  phase  of  the  carrier  is  shifted  hy 
180°  (figure  7). 


wn/wnAAA/vn/un*- 

Figure  7.  Biphase  Modulated  Carrier. 

The  basic  lengths,  periods,  and  frequencies  describing  the  carrier,  pulses,  and  code, 
together  with  the  other  parameters  used,  are  given  in  Table  1. 

Table  I. 

Length  Period  Frequency 


Carrier 


Pulse- 


String 


1.903429  x 10-*  m 

6.349841  x 10" 10  sec 

1.575  x lCT*  Hz 

/~  20  cm\ 

8 in.  > 

(~  6 nanoseconds) 

<1  575  GHz) 

2.997925  x I0‘  m 

10'7  sec 

1.0  x I07  bit/sec 

/ ~ 30  m \ 

/ ^ — 0.1  microsecond  \ 

(10  megabits/sec) 

\ ~ 33  yd/ 

V ~ 1 00  nanoseconds/ 

3.0669  x I04  m 

1.023  x I0"4  sec 

9.7751 7 x I03  sec'1 

/ — 31  km\ 

(~  0. 1 milliseconds) 

(~  10,000/sec) 

\~  19  mi  / 

Signal  Velocity  = C = 2.997925  x 10*  m s'1 
Pulse  Length  21  1 57.5  Wavelengths 
Code  Length  = 1 ,023  Pulses 

The  various  parameters  given  are  important  in  determining  the  ranges  of  path  length 
difference  over  which  the  different  mechanisms  to  be  described  become  operative. 

As  an  example,  consider  the  case  when  the  strongest  and  most  direct  component 
of  the  incident  signal  hits  an  obstacle,  such  as  a trunk  or  branch  of  a tree.  The  obstacle 
may  diffract  the  signal  into  a diffraction  pattern  having  a characteristic  spatial  intensity 
distribution.  At  certain  locations  of  the  pattern,  the  intensity  may  be  very  small.  A 
receiver  that  is  locked  to  this  signal  may  lose  it  if  the  receiver  antenna  is  moved  to  a 
low  intensity  point  of  the  diffraction  pattern.  If  there  is  another  sufficiently  strong 


hut  delayed  component,  the  receiver  may  acquire  it  at  this  point  and  the  resulting 
range  measurement  would  produce  a faulty  range 

The  intensity  variations  in  the  diffraction  pattern  can  he  explained  by  the 
processes  of  constructive  and  destructive  interference.  I his  type  of  interference  is 
observed  when  coherent  signals  or  waves  are  superimposed  However,  incoherent  sig- 
nals or  waves  do  not  produce  observable  interference  patterns  Strictly  speaking, 
destructive  interference  doesn’t  exist  for  pseudo-noise  biphase-modulated  signals 
because  the  phase  of  the  carrier  changes  almost  randomly  by  1X0°  in  accordance  with 
the  modulation.  If  the  path  difference  of  two  components  is  larger  than  the  modula- 
tion pulse  length,  the  interference  at  a given  point  changes  on  the  average  with  each 
clock  pulse.  Because  of  the  very  high  clock  rate,  no  destructive  interference  can  be 
observed,  l or  path  differences  less  than  the  modulation  pulse  length,  destructive 
interference  may  be  observed.  If  (i  is  the  path  difference  and  L the  pulse  length  and 
0 < C < L,  the  observable  portion  of  the  interference  is  L - 8/L,  and  the  nonobservable 
portion  is  li/L. 

However,  if  two  signals  were  present,  but  separated  in  path  length  by  about  a 
pulse  length  (30  meters)  or  more,  the  tendency  toward  carrier  cancellation  would 
depend  on  the  autocorrelation  function  for  the  code.  The  actual  codes  used  are  Gold 
Codes  possessing  pseudorandom  characteristics;  therefore,  beyond  a relative  path 
difference  of  about  a pulse  length  (30  m),  the  signals  tend  to  interfere  constructively 
about  as  often  as  destructively,  and  the  mechanism  of  carrier  cancellation  is  inoperative. 
In  the  past,  the  combination  of  (1)  normal  operation  in  a nonscattering  environment, 
and  (2)  over-emphasis  of  the  dependence  of  the  signal  correlation  on  the  code  correla- 
tion properties  led  to  the  neglect  of  possible  carrier  interference  effects  when  operating 
in  a scattering  or  multipath  environment.  Consequently,  the  importance  of  the 
phenomenon  of  carrier  cancellation  was  underrated. 

The  actual  codes  used  in  the  simulations  were  shorter  than  the  codes  normally 
used  in  the  field,  but  they  possessed  similar  characteristics.  The  code  generation 
system  that  was  simulated  consisted  of  a pair  of  shift  registers,  connected  as  shown  in 
figure  8.  Bach  shift  register  consists  of  10  consecutive  binary  storage  elements  called 
stages  and  a feedback  network.  The  binary  values  (c.g.  two  discrete  voltages)  that  can 
be  stored  in  each  stage  are  for  the  purpose  of  this  treatment  defined  as  value  zero  (0) 
and  value  one  (+1 ).  All  stages  of  the  two  shift  registers  in  figure  6 show  the  value  (+1 ). 
The  feedback  network  feeds  the  values  of  some  designated  stages  into  a modulo  two 
adder.  Depending  on  the  input  values  of  the  adder,  the  output  of  the  adder  will  assume 
the  values  0 or  +1.  The  output  of  the  adder  is  applied  to  the  input  of  the  shift  register. 
The  outputs  of  both  shift  registers  are  fed  into  a modulo  two  adder  that  generates  the 
Gold  Code  modulation  s(t).  If  the  values  at  the  output  of  both  shift  registers  are 
equal,  a value  of  (+1)  is  assigned  to  the  code;  if  they  are  unequal,  a value  of  (-1)  is 
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PULSE 


SHIFT  RESISTED  * \ 
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Figure  8.  Gold  Code  Generation. 


\ assigned  to  the  code.  The  shift  registers  are  operated  by  applying  a clock  pulse  to  each 

stage.  The  application  of  a clock  pulse  causes  the  value  of  any  stage  to  be  shifted  to 
the  next  stage  on  the  right  and  the  output  of  the  feedback  network  to  be  shifted  into 
the  first  stage  of  the  shift  register.  To  initialize  the  registers,  each  element  of  the 
registers  was  set  to  +1.  Provision  was  made  in  the  computer  program  for  representing 
different  feedback  connections  (only  certain  ones  would  be  capable  of  generating  a 
maximum  length  code).  Also  one  of  the  registers  could  be  advanced  independently  of 
the  other  to  change  the  code  and  to  determine  the  effects  of  such  changes.  An 
example  of  one  of  the  Gold  Codes  actually  used  is  given  in  figure  9.  The  autocorrela- 
tion function  for  this  code  is  given  in  figure  10.  Its  principal  correlation  value  is  1 ,023 
and  is  determined  by  the  integration  range  that,  in  this  case,  is  equal  to  the  code  length. 
The  other  values  assumed  by  the  correlation  function  are  63,  -65,  and  -I.  It  would 
seem  desirable  to  use  a code  that  has  correlation  values  of  - I as  the  only  value  occur- 
i ring  near  the  principal  correlation  value  of  1 ,023.  Some  simulations  have  shown  that  it 

the  code  is  improved  in  this  respect  by  advancing  one  shift  register  relative  to  the 
other,  the  ability  to  correlate  can  be  improved. 

MATHKMATICAL  ANALYSIS 

Correlation  of  the  Received  Multicomponent  Signal  with  a Reference  Signal. 

To  perform  the  simulations,  it  was  necessary  to  derive  the  computational  equations 
and  to  formulate  a method  to  reduce  the  computations  to  a tractable  size.  Also,  the 
resulting  equations  should  lend  themselves  to  a convenient  interpretation  of  the 
problem.  The  signal  was  assumed  to  be  repeated  periodically  without  pause.  The 
equations  were  presented  in  the  time  domain  rather  than  in  the  spatial  domain. 


8 


Figure  9 Example  of  a Gold  Code 


i«Uf  1*1  ION  r W>H.  1 t J* 


• 1 

41 

-1 

-i 

•43 

•1 

61 

*1 

•1 

• 1 

• i 

•3 

-1 

-1 

• 3 

•3 

-1 

•3 

•1 

63 

•46 

-63 

- 1 

M 

-t 

41 

•3 

• 3 

-•3 

•1 

•3 

-3 

4J 

-1 

-3 

-3 

•3 

•1 

%1 

-I 

-1 

-3 

-t 

44 

41 

-1 

•1 

-4* 

-1 

-3 

•3 

-43 

-1 

-4 

-3 

•1 

-1 

41 

•3 

• 1 

• | 

•3 

-4 

-1 

-1 

• *. 

-3 

•3 

•3 

•63 

-3 

•4 

-1 

• 3 

• 1 

-3 

•3 

• 3 

6 1 

-4 

-1 

•1 

•1 

• 4* 

41 

61 

-1 

-1 

-1 

-3 

-3 

•3 

44 

-3 

•4 

-» 

-1 

-l 

-3 

-3 

• | 

•3 

-4 

•1 

o3 

-3 

*1 

-43 

-63 

•4 

-1 

-1 

-3 

41 

•3 

-3 

-4 

-l 

-l 

-3 

•3 

-3 

-l 

-4 

•1 

•X 

* 3 

43 

61 

* 9 

44 

-1 

•4 

-1 

•X 

-*t» 

-3 

-3 

• | 

64 

-4 

-1 

•X 

-3 

-3 

-3 

41 

-43 

-3 

-4 

-1 

-l 

-3 

-3 

-3 

-63 

-1 

61 

61 

-I 

-33 

-3 

-43 

-63 

44 

61 

31 

-1 

-I 

-1 

-3 

-3 

-1 

-4 

-I 

-l 

-43 

-3 

-1 

- t 

-43 

-63 

-4 

-1 

-43 

-3 

41 

•3 

-63 

44 

• | 

-63 

•4 

-1 

-1 

• 1 

41 

-3 

-43 

-43 

-4 

-1 

-69 

43 

-3 

-63 

-43 

-1 

-4 

-1 

-I 

-3 

-3 

-3 

64 

-4 

-4 

-1 

-l 

-3 

43 

-3 

44 

-1 

-4 

61 

-I 

-3 

41 

•3 

4] 

-l 

-4 

ft* 

41 

4 J 

-1 

4J 

-43 

-1 

-4 

-I 

-1 

41 

-3 

-3 

41 

-4 

-1 

-l 

-3 

-43 

-1 

-I 

-4 

-3 

-43 

-3 

-l 

-1 

44 

-1 

-4 

-I 

-1 

43 

-3 

-3 

-63 

-43 

-1 

-4 

-1 

-1 

6! 

61 

-» 

-43 

-I 

-1 

-4 

43 

-1 

-1 

•1 

44 

-1 

-4 

-1 

-1 

-1 

-l 

65 

-1 

-63 

-l 

-l 

-1 

-3 

-1 

44 

44 

-l 

-1 

-1 

-43 

-3 

-43 

-l 

-4 

-4 

64 

4? 

-l 

-1 

-3 

44 

4 J 

64 

bS 

-3 

-3 

31 

-4 

44 

-» 

-43 

-3 

-3 

-1 

-I 

-4 

-1 

-1 

-I 

-43 

41 

-1 

44 

-43 

-1 

-1 

-l 

-3 

-1 

-1 

-I 

-I 

-4 

-43 

-I 

43 

-3 

-1 

-1 

44 

-4 

44 

*3 

43 

43 

-1 

-1 

-4 

64 

-1 

-I 

-3 

-3 

-1 

-43 

-4 

-63 

-1 

-1 

-3 

-3 

• 1 

61 

•4 

-4 

-l 

-43 

-3 

-t 

-3 

44 

-4 

-4 

-l 

-1 

-3 

63 

-3 

-4 

41 

-1 

-l 

-3 

-3 

-33 

41 

41 

-l 

6 3 

-l 

-3 

-63 

“3 

-63 

44 

-l 

-1 

-l 

-3 

-3 

33 

• 3 

-4 

-64 

-1 

-l 

-3 

-3 

-1 

-1 

•63 

-1 

-i 

-1 

-1 

43 

-3 

-1 

-1 

-1 

44 

-1 

b3 

-l 

-3 

-3 

-63 

31 

41 

64 

-1 

-1 

-l 

-1 

-3 

- 1 

64 

64 

-4 

-4 

-1 

-1 

-3 

61 

*3 

- 3 

61 

-4 

-I 

-1 

-1 

-3 

-1 

-3 

61 

-4 

-l 

-1 

-I 

-3 

-3 

-63 

• I 

-43 

41 

-63 

-1 

-l 

-43 

-63 

-3 

• 1 

44 

41 

-4 

bl 

-1 

-43 

-3 

61 

-63 

44 

-4 

•hi 

•3 

-l 

S 

-63 

-3 

-1 

-43 

-1 

-1 

-1 

• ] 

-1 

-1 

-4 

-1 

-1 

43 

4 J 

-1 

6 9 

-43 

•43 

-4 

-69 

-1 

43 

4) 

-3 

-1 

-1 

-43 

-4 

-4 

-1 

-1 

-1 

-43 

-1 

4) 

-4 

-4 

-1 

-1 

43 

-1 

*1 

-43 

-4 

-1 

-1 

-1 

-1 

61 

44 

44 

41 

4 J 

-4 

-1 

-3 

-3 

-i 

-1 

- 1 

-4 

-4 

-1 

•l 

-I 

-I 

-1 

44 

-4 

-4 

-l 

-1 

-43 

-i 

•I 

-43 

44 

• 1 

44 

-4 

-1 

-1 

-3 

-3 

-3 

-4 

-4 

43 

-1 

-1 

41 

•1 

-4 

-1 

-l 

•l 

41 

- 1 

•l 

f 1 

-63 

-43 

•4 

-4 

-» 

-1 

41 

-I 

-3 

• | 

•4 

64 

-I 

-1 

-3 

•43 

•3 

- 1 

•4 

-4 

-1 

•l 

-1 

•« 

-I 

44 

-4 

-4 

-1 

-1 

-l 

-t 

-33 

- 1 

-1 

-63 

•1 

-1 

43 

-1 

-3 

61 

-4 

-4 

-I 

-l 

•l 

-1 

-1 

•4 

64 

-I 

-1 

-1 

-43 

•3 

61 

-4 

64 

-1 

-43 

-43 

41 

-1 

•I 

-1 

-4 

-I 

• 1 

-1 

• I 

-1 

*1 

•4 

-4 

-1 

4 J 

43 

1429 

-|  *t  •»  -»  • 3 -43 

Figure  10.  Autocorrelation  Function  for  a Gold  Code. 
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The  signal  produced  hy  the  transmitter  is  given  by 


S(t)  = A s(t)  cos  cot 


(I) 


\ 


i 


i 


I 


where  A is  the  amplitude.  u>  is  the  carrier  frequency,  t is  the  time  and  stt)  is  the 
modulation  representing  the  code  (s(t)  = 1 1).  The  signal  at  the  receiver  S(  (t)  is  the 
superposition  of  a number  M of  component  waves  S(  (t);  each  being  a replica  of  the 
transmitted  signal  and  having  an  amplitude  and  a time  delay  r . 


S,  (t) 


S,  (t) 


E 

j=i 


rvi 

E 


S.  (t) 


)*> 


Aj  s(  t - r( ) cos  co(  t - t ) 


(2) 

(3) 


The  delays  r i must  account  for  (I)  increased  travel  distance,  (2)  altered  velocity  of 
propagation,  and  (3)  phase  changes  due  to  scattering.  For  purposes  of  computation, 
they  are  expressed  in  terms  of  equivalent  ranges  (optical  paths)  r|t 


fj  = c r.  (4) 

where  c is  the  velocity  of  propagation  in  a vacuum. 


At  the  receiver,  the  modulated  signal  is  separated  from  the  carrier.  This  is  repre- 
sented mathematically  by  multiplying  the  incoming  signal  Sr(t)  with  an  internally 
generated  signal  SR  (t)  = Aq  cos  co(t  - ro).  The  demodulated  signal  Sm(t)  is 


M 

Sm(t>  = S,(t)  SR(t)  = Ao  A.  s(t  - r ) cos  oj( t - r ) cos  co( t - rQ ) (5) 


j*l 


M 

sm  (t)  = ^ Ao  A;  s( t - t.  ) j Zi  cos  cu(2t  - Tj  - ru ) + Vi  cos  - rQ ) | (6) 


1*1 


In  practice,  demodulation  may  be  accomplished  by  a somewhat  different  procedure, 
but  the  result  will  correspond  ideally  with  the  mathematical  results  given  here.  The 
term  in  (6)  with  the  factor  '/2Cosco(2t  - r - ro)  = Vi  cos(2wt  - <p)  represents  the  high 
frequency  component  with  the  angular  frequency  2w.  This  component  can  be  filtered 
out  by  a low  pass  filter.  The  other  component  of  (6)  represents  the  detected  signal. 


1 1 


u. 


(7) 


I 


or 


with 


Sa  = 


M 

s 


J"1 


COS  U)(T)  - To)  S(t  - Tj) 


M 


sd  = SB's"'T,) 

j'l 

Ao  A, 

_j=  co_,.j  ,o. 


(8) 


Equations  (3)  and  (7)  are  suitable  to  demonstrate  that  carrier  cancellation  is  a property 
of  the  scattered  signal  and  cannot  be  influenced  by  receiver  manipulations.  Consider 
the  case  that  the  signal  Sf(t)  at  the  receiver  consists  of  the  three  components: 

S,(t)  = A,  stt-T,)coscu(t-T,) 


S2(t)  = A2  s(t  - r2 ) cos  ou(t  - r2 ) 


S3(t)  = A3  s(t-r3)coscj(t-r3) 

With  t2  = t,  + — r , and  f = carrier  frequency,  S|  (t)  and  S2(t)  are  the  components  of 

the  direct  signal,  and  S3(t)  is  the  delayed  component.  Because  the  pulse  duration  T of 
the  modulation  is  large  compared  to  — , S(t  - r( ) is  equal  to  S(t  - t2  ),  except  in  the 

vicinity  where  S,(t)  or  S2(t)  change  signs.  For  the  cosine  terms  of  S,(t)  and  S2(t), 
it  is 

cos  u>(t  - t2  ) = cos  (cut  - cut,  - oj/20  = cos  I co< t - r , ) - rr | = - cos  a)(t  - t,  ). 

Therefore,  S,(t)  and  S2(t)  cancel  out  except  in  the  vicinity  where  S,  (t),  or  S2(t) 
change  signs.  The  remaining  component  S3(t),  on  which  the  receiver  may  lock,  will 
cause  a range  error  that  is  proportional  to  t3  - rt . 

In  the  detected  signal  of  equation  (7),  only  AQ  and  tq  are  determined  by  the 
receiver.  Normally  ro  is  adjusted  by  the  receiver  such  that  cos  cu(r2  - rQ)  = 1. 

The  detected  signal  is  to  be  correlated  with  a reference  signal  that  has  the  form 


I 


12 


K(t  - 6>  = A(  s< t - 6)  (V) 

where  6 is  the  phase  of  the  modulation.  It  should  be  emphasized  that  the  phase  of  the 
reference  carrier  ro  and  the  phase  of  the  reference  modulation  6 must  be  retained  as 
independently  variable  quantities,  if  a realistic  simulation  is  to  be  obtained  The  cross- 
correlation  function  06)  between  the  received  and  the  reference  signal  is  to  he 
evaluated. 

06)  = lim  f R (t  - 6)  S.  (t)  dt  (10) 

L -*«>  2L  J a 

■ i 


The  value  of  6 that  maximizes  it  will  give  the  time  delay  and,  consequently,  the  range. 

The  signal  s(t)  is  composed  of  a series  of  pulses,  each  having  duration  T,  that 
make  up  a coded  string  of  N pulses  having  duration  NT.  By  introducing  the  function 
a (x)  given  on  the  interval  (0,  N)  by 


a (x) 


I;  0 <x<  I 
* 0;  1 < x < N 

a single  pulse  assumes  the  following  form: 

1 ; 0 < t<  T 
: < NT 


(II) 


v ( I;  0<  t« 

°(r>  !0;t<, 


(12) 


with  — = x. 
T 


Outside  of  this  fundamental  interval  0 < x < 1 , the  function  will  be  defined  as 
periodic  of  period  N: 

(13) 

(14) 

where  k assumes  positive  and  negative  integral  values.  The  function  o (x)  is  shown  in 
figure  1 1. 


a (x  + k N)  = a (x) 


/ t + k NT\  / t \ 
0 ( T / ° (f ) 


l 
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Figure  11.  The  Function  o(x|. 
The  coded  string  of  pulses  now  assumes  the  form 


s(t)  = a(,  j<Y<j+,’aja±l  f lf>) 

the  actual  values  of  a(’s  being  given  by  the  code.  The  periodicity  of  o accounts  for  the 
fa  t that  s(t)  is  transmitted  periodically: 

s(t  + kNT)  = s(t),  k = O.i  l,±  2 < 1 7) 


The  defected  signal  and  the  reference  signal  may  now  be  expressed  in  more  explicit 
terms  in  the  forms  given  by  equations  ( 18)  and  ( I1)),  respectively. 


R(t-S> 


(18) 


<ll>> 


Because  of  the  periodicity  of  its  integrand,  the  expression  for  the  correlation  function 
given  by  equation  (10)  may  be  replaced  by  one  having  an  integral  over  only  one 
period ; 

NT 

C(6)  = ~ f R(t  - 6)  Sd(t)  dt  (20) 

O 


14 


M N N 


i=  1 )=  1 k = 1 


By  making  the  substitution 


t'  = t - 6 - kT, 


equation  (21 ) can  be  expressed  in  the  following  form: 


A 

C<8)  = — 
NT 


M N N 


NT-6  -kT 


Bi  aj  ak 


i*  1 j=l  k = 1 -6-kT 


f / t'  + 6 + kT  - r - jT  \ /t. 

J ° T-1—  “(f) d'  (-3* 

LT  ' ' 


Again  because  of  periodicity  of  the  integrand,  the  integral  may  be  taken  over  any 
complete  period: 


M N N 


•“  * i->  k = i 

Using  equation  ( I 2)  to  define  a in  the  second  factor  of  the  integrand,  the  form  of  the 
integral  may  be  further  simplified: 


t'  + 6 - Tj  - (j-k ) T\  /, 


) 04, 


M N N 


/ /I'  + S-r.-U-klT1,  , 

“iVi  / f / d' 


i-l  j-i  k = l 


If  vi  is  taken  to  be  the  largest  integer  such  that 


"*  < T 


and  Xj  is  taken  to  be  the  fractional  remainder 


-&  + T, 

Xi  T ‘ 


0 < Xj  < I 
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O',  + xj)T  = -6  + r| 


the  correlation  function  assumes  the  form 

M N N ^ 


('(6) 


NT 


I )*>  k- 


i r 


or  letting  t”  = — 


06) 


M N N 

SY 2 B'a'a‘T O-k-x, )<!.■' 


i=  1 j=  I k ■ I 


( 


(30) 


(31) 


From  the  actual  physical  situation,  it  is  known  that  signal  components  lor  which  the 
inequality 


0<r,«NT  (32) 

does  not  hold  will  give  only  an  insignificant  contribution  to  the  received  signal,  since 
otherwise  the  signal  will  have  traveled  over  such  a long  path  (with  possible  multiple 
scattering)  that  it  will  be  severely  attenuated.  Thus,  for  correlation,  the  only  realistic 
phases  for  the  reference  signal  correspond  to 

0<  6 « NT  (33) 


or 


I y-  + I < < N (34) 

In  evaluating  the  integral,  values  of  vi  need  only  be  taken  in  the  following  range: 

I I < < N (35) 

As  an  example,  more  than  500-meter  foliage  penetration  will  completely  destroy  the 
incident  wave;  hence,  I vi  I < 17.  Within  this  range  of  values  of  v.,  the  integral  in  the 
expression  forC(6)  is  readily  evaluated  (figures  12  and  13). 

By  introducing  a periodic  analog  of  the  Kronecker  6 defined  by 


16 


V 


il ; i - j = nN,  n = 0,  ± I , ± 2. . . 
0;  i - j is  not  a multiple  of  N 


(3ft| 


the  value  of  the  integral  may  he  written  in  the  form  given  below 
t 


/ 


o (t”  - v,  + j - k - X,)  dt”  = X,  + <l-\>  fi^.^k.0 


(37) 


The  correlation  function  C'(6)  given  in  equation  (31)  may  now  be  expressed  in  the 
form 


m n n 


A 

((6)  =»  v 
N 


B.a,ak  [Xi6j.vj+k+t  +(,-Xi)6,.,j+k]  (38> 


i*  1 j=  I 


and  by  periodically  extending  the  range  of  subscripts  of  af  such  that 

aj+N  = aj’ 

summation  over  j may  be  performed: 

M N 

c(6)=  h [Xj +(|-w*ki 

1-1  k=  I 

The  autocorrelation  function  for  the  code  is  defined  by 

N 

7«  = ^ ak  a«*k 

k - 1 

and  may  be  used  in  the  expression  for  C( 5): 

M 

081 ' £ 23  B‘[vx'‘V.'V] 

i-  I 


(39) 


(40) 


(41) 


(42) 


If  the  7#’s  are  calculated  once  and  stored  in  a computer,  they  may  be  used  repeatedly 
in  the  computation  of  C’(6). 
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Adjustment  of  the  Receiver’s  Carrier  Phase  and  Liquations  for  Computing  the 
C orrelation  Function.  By  using  equation  (8)  and  hy  introducing  the  coefficients 


*1- 


•V\A 

2N 


x (7 


V 


<4.0 


an  alternative  expression  for  the  correlation  function  may  be  written  in  a form  where 
the  carrier  phases  of  the  component  waves  is  emphasized  ' 


06)  = 


0.  cos  to  (r( 


(44) 


If  the  carrier  phase  to  of  the  receiver  generated  signal  is  considered  to  be  variable,  but 
all  of  the  other  parameters  in  equation  (44)  are  fixed,  then  the  correlation  function 
06)  may  be  written  as 

C(8)  = C cos  to  (f-  to  ) (45) 


since  a sum  of  sinusoidal  functions  of  a single  frequency  is  again  a sinusoidal  function 
of  the  same  frequency.  At  the  receiver,  6 is  shifted  to  find  the  maximum  value  of 
06),  but  as  6 is  shifted  ro  also  follows  in  such  a way  that  06)  is  maximized  in  equa- 
tion (45)  so  that  06)  = C and  td  = ¥ + 2mr  (where  n = integer).  Again,  this  is  an 
idealization,  and  the  design  and  construction  of  actual  receivers  may  correspond  only 
in  this  ideal.  In  actuality,  6 is  shifted  by  the  voltage  controlled  oscillator  (VCO)  that 
drives  the  code  generator,  and  to  is  adjusted  by  the  voltage  controlled  oscillator  (VCO) 
of  the  carrier  phase  lock  loop.  After  signal  acquisition,  code  tracking  loop  and  phase 
lock  loop  may  operate  coherently.  From  equations  (44)  and  (45), 

M 

0,  cos  to  (t;  - rH)  (46) 

i*  1 

and  the  functional  form  of  the  equation  must  be  preserved  for  any  translation  of  the 
to  axis: 

M 

sin  to  (ri -ro)  (47) 

(•  l 

Also,  for  the  same  reason,  equations  (46)  and  (47)  must  hold  for  the  particular  choice 

To  =0’ 


C sin  to  (¥-  t ) = 


E 


C cos  to  (7-  to  ) = 


E 


C sin  u)  t = 


■ 


pt  sin  oj  r, 


M 

C cos  u>  f= 

i*  i 


cos  u)  r, 


and  from  these  relationships,  C is  readily  determined: 


(48) 


(49) 


CJ  = (rsincaT)3  +(CrcoswT)J 


(50) 


r = 


(51) 


Equation  (51)  was  used  in  the  calculation  of  the  correlation  function  ('  by  computer. 
If  desired,  the  phase  Tcan  be  found  from  equations  (48)  and  (49)  if  care  is  taken  in 
removing  phase  ambiguities  that  can  arise  in  calculating  inverse  trigonometric  functions. 


The  Maximum  Value  of  the  Correlation  Function.  Although  equation  (51)  will 
suffice  for  computing  ^(6),  computing  this  function  at  enough  points  to  insure  finding 
its  maximum  can  still  present  a formidable  problem.  However,  by  introducing  a com- 
plex function  C(6)  and  by  dividing  the  6-axis  into  certain  intervals,  an  interpretation 
can  be  made  that  will  permit  P(6)  to  be  computed  only  at  a finite  and  reasonable 
(although  large)  number  of  predetermined  points.  IfC'(6)  is  given  by 


C (6) 


k » 1 


(52) 


then  by  using  de  Moivre’s  relation  for  the  complex  exponentials. 
M M 

C(6)  = 0k  cos  u>  rk  + i ^ Pk  sin  cj  rk 


(53) 


k - I 


k * 1 


it  is  evident  that  C( 6)  is  the  modulus  or  distance  from  the  origin  to  the  point  C(6) 
in  the  complex  plane: 


C(6)  = IC(6)I  = |C(6)C*(6)|* 


(54) 
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I 


11  the  quantities  xj  are  taken  to  be  a permutation  of  the  quantities  x(  such  that 

x*  < x'  < < x , (55) 

12  n 

then  the  particular  values  of  6 denoted  by 

6m  M n = (x'  + n)T  (56) 

m ▼ M n m 

will  have  the  following  property. 

«„<«».  i 1571 


Furthermore,  examination  of  equation  (43)  and  (53)  reveals  that  between  successive 
values  of  6k  , C(6)  is  a linear  function  of  S and  will  move  in  the  complex  plane  along 
the  straight  line  joining  the  two  points  C(8k  ) and  C(5k+1 ) (tigure  14), 


5-\  r- 

C(6)  = C(6k  ) + — |C(6k  + , 

6k  + i-6k  L 


-C(6k)];6k<6<6k+1 


(58) 


c - plane 
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Within  each  interval  defined  in  this  way,  05)  can  have  at  most  one  minimum  value, 
and  its  maximum  value  must  occur  at  one  of  the  endpoints.  It  can  easily  be  shown 
that  on  any  one  of  the  intervals,  05)  is  represented  by  a conic  section.  The  endpoints 
of  the  intervals  are  the  points  that  have  been  calculated  in  the  computer  program. 
Thus,  for  a code  of  length  N and  for  M components,  the  correlation  function  C(8) 
needs  to  be  computed  only  at  NM  points  to  insure  finding  its  maximum. 

The  Incoherent  Limiting  Case.  In  the  case  where  a very  large  number  of  com- 
ponent waves  are  present  and  all  have  about  equal  amplitude  and  randomized  phases, 
the  correlation  function  can  be  shown  to  be  independent  of  the  carrier  phase,  equa- 
tion (51)  may  be  rewritten  in  the  form 


M m 


i=i  j=i 


or,  by  using  a familiar  trigonometric  identity,  in  the  following  form: 


C2 


p.  p.  cos  OJ  (Tj  - Tj) 


(60) 


In  equation  (60),  terms  for  which  i is  different  from  j will  occur  in  pairs  with  the  roles 
of  i and  j interchanged;  and  since  the  value  of  the  cosine  is  independent  of  the  sign  of 
its  argument,  the  second  summation  can  be  replaced  by  an  expression  involving  fewer 
terms. 


i-l 


C2  = 


i»  1 


' I j=l 


(61) 


(62) 


In  the  case  of  a large  number  of  components  with  randomized  relative  phases,  the 
second  term  will  tend  to  be  small  relative  to  the  first  since  the  positive  terms  can  be 
cancelled  by  negative  ones;  whereas,  in  the  first  term  all  contributions  are  positive. 
In  this  incoherent  limiting  case  or  random  phase  approximation,  the  correlation  func- 
tion becomes  independent  of  the  carrier  phases  of  the  components 


22 


06)  = 


(63) 


and  this  corresponds  with  the  presence  of  many  multipath  signals  of  about  equal  ampli- 
tude whose  path  lengths  are  randomized  relative  to  a wavelength. 

DISCUSSION  OF  COMPUTED  EXAMPLES 

By  using  the  formulas  and  methods  previously  outlined,  a computer  program  was 
written  and  the  results  were  computed  for  examples  that  seemed  to  conform 
reasonably  with  the  experimental  situation  and  that  would  demonstrate  the 
phenomenon  of  carrier  cancellation.  Examples  of  correlation  functions  are  given  in 
figures  15  and  16.  Both  examples  show  results  for  seven  components  with  the  assigned 
amplitudes  and  path  lengths  shown  below: 

Range  (meters)  Amplitude 


200.0 

m 

1.00 

201.0  + 6 

in 

0.95 

230.07442 

m 

0.20 

260.05367 

m 

0.20 

290.03292 

m 

0.20 

320.01217 

m 

0.20 

349.99142 

in 

0.20 

For  the  correlation  function  in  figure  15,  the  path  increment  6 corresponded  to  270° 
in  carrier  wavelength,  and  in  figure  16,  6 corresponded  to  90°. 

The  data  are  plotted  again  in  figures  17  and  18  for  comparison  purposes.  Here, 
suggestive  connecting  curves  have  been  added,  although  their  exact  form  is  not  known. 
The  change  at  around  200  meters  is  a clear  example  of  carrier  cancellation.  In  this 
example,  either  one  of  the  first  two  components  contain  more  than  eight  times  the 
combined  energy  of  the  remaining  five  components;  yet,  in  the  case  of  carrier 
cancellation,  the  latter  five  components  all  show  correlation  values  at  least  twice  as 
high  as  that  of  the  two  larger  components.  This  all  occurs  because  the  path  length  of 
one  of  the  strong  components  is  changed  by  half  a wavelength,  or  about  1 0 centimeters. 
It  should  be  noted  that  the  five  minor  components  shown  appear  to  be  almost 
completely  unaltered  by  the  large  changes  associated  with  the  two  principal  com- 
ponents. The  reason  that  they  remained  unaltered  is  that  the  minor  components  are 
separated  from  the  principal  components  in  path  length  by  a pulse  length  or  more. 
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Figure  15  First  Correlation  Function. 


NORMALIZED 

CORRELATION 

FUNCTION 


Figure  16.  Second  Correlation  Function 
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Figure  17.  First  Correlation  Function  with  Connecting  Curve. 
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Figure  18.  Second  Correlation  Function  with  Connecting  Curve 


The  change  in  correlation  peaks  can  be  followed  as  the  phase  of  one  of  the  com- 
ponents is  changed,  l or  example,  the  curves  in  figure  19  are  plotted  lor  three  com- 
ponent waves  corresponding  with  the  following  ranges  and  amplitudes: 


Range  (meters)  Amplitude 


300.0 

m 

1.0 

300.0  + 6 

m 

1.0 

400.0 

m 

0.5 

In  this  case,  the  abscissa  gives  the  increase  in  the  path  length  6 associated  with  com- 
ponent 2.  The  correlation  peaks  that  are  followed  correspond  to  300  and  400  meters, 
respectively.  The  complete  cancellation  in  the  300-meter  correlation  peak  occurs 
because  the  first  two  components  have  exactly  the  same  amplitude  and  almost  no 
range  separation. 


IV.  CONCLUSIONS 


It  is  concluded  that 

1.  The  proposed  mechanism  of  carrier  cancellation  is  a plausible  explanation 
for  anomalous  range  measurements  obtained  in  a scattering  environment  using  a code 
modulated  signal. 

2.  Carrier  cancellation  is  a property  of  the  scattered  signal  only  and  cannot  be 
influenced  by  receiver  signal  processing. 

3.  Carrier  cancellation  may  occur  when  the  path  difference  of  the  interfering 
signals  is  less  than  the  length  of  the  modulation  pulse. 

4.  In  the  incoherent  limiting  case  where  a large  number  of  component  waves 
are  present,  the  correlation  function  becomes  independent  of  the  carrier  phases  of  the 
components.  Then,  a receiver  will  lock  onto  the  signal  only  when  the  correlation 
function  has  distinctive  peaks. 

5.  The  probability  of  large  ranging  errors  due  to  carrier  cancellation  is  inversely 
related  to  the  modulation  bit  rate. 
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APPENDIX 


THE  COMPUTER  PROGRAM 


The  computations  were  performed  on  a ( DC  6600  computer  using  a program 
written  in  FORTRAN.  The  program  consisted  of  a main  program  and  three  sub- 
routines. Subroutine  CODEIN  is  used  to  generate  the  Gold  Code  and  to  store  it  in  a 
COMMON  area.  Its  argument,  NADVA,  specifies  the  relative  advance  of  the  two  shift 
registers.  The  amount  of  relative  advance  may  be  changed  without  altering  the  con- 
figuration of  shift  register  connections  by  using  ENTRY  SHFTADV.  The  autocorrela- 
tion function  for  the  code  is  computed  by  calling  subroutine  CORVAL,  which  in  turn 
computes  each  value  of  the  autocorrelation  function  by  calling  subroutine  CODECOR. 
The  main  program  computes  the  correlation  function  of  the  multicomponent  signal 
with  the  receiver  generated  signal.  This  is  computed  at  all  of  the  predetermined  points 
(mentioned  previously).  The  correlation  function  is  printed  out  in  detail  for  the 
principal  range  (or  range  differences)  of  interest.  For  the  rest  of  the  correlation  func- 
tion, the  largest  values  and  their  corresponding  ranges  are  stored,  and  only  these  are 
printed  out.  A list  of  input  data  with  the  corresponding  card  formats  is  given  below 
followed  by  a listing  of  the  actual  computer  programs. 
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INPUT  DATA 


Field 

Card 

Specification 

Column 

CARD  1 NCC 

Number  of  Computed  Simulations 
to  be  Run 

15 

1-5 

FREQ 

Frequency  of  Carrier  (Hertz) 

El  5 

6-20 

PULE 

Pulse  Frequency  ( 1 .023  E+07 

El  5 

21-35 

(bit/sec)) 

VEL 

Velocity  of  Propagation  (Meters/ 

E 1 5 

36-50 

Second) 

CARD  2 Contains  shift  register  configuration  for  the  two 

1011 

1-10 

shift  registers  used  in  generating  the  code,  e.g. 

1011 

21-30 

00 1 000000 1 , 0 1 1 00 1 0 1 1 1 specifies  connections 
for  feedback  of  the  sum  (mod  2)  from  the  3d  and 
1 0th  registers  of  the  first  shift  register  and  from 
the  2d,  3d,  6th,  8th,  9th,  and  10th  register  in  the 
second. 


NCC  Card  Sets  for  each  Complete  Simulation 


CARD  1 

NPW  - Number  of  Partial  Waves 

15 

1-5 

NC’RV  — Number  of  Correlation  Values  to  be 
Computed 

15 

11-15 

IADVA  — Relative  Advance  of  th“  first  Shift 
register  with  respect  to  the  second 

NPW  Card  Sets  for  each  Partial  Wave  (One  Card  per  Set) 

15 

21-25 

CARD  1 

PL  - Path  Length  for  the  Partial  Wave 

(meters) 

El  5 

1-15 

AMP  - Amplitude  of  the  Partial  Wave 

(arbitrary  units) 

F5 

21-25 

NPI  - Number  of  Phase  Increments 

15 

31-35 

TPI  - Total  Phase  Increment  or  Range  of 

Phase  Angles  to  be  Covered  (degrees) 

F5 

41-45 
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t PROGRAM  C0UC3RIINPUT, OUTPUT, TAPE5*lNPUT,TAPFb*0UIPuTl 

COMMON  ICI102  3) . IC»TL(1023) , I A lldl  .13(101 

DIMENSION  °L(10I.AMP(10).NPI (10l,TPI(10).CHl(10) , CHlOl 10  I . 1 / 'll  1 t 
•.BESMI101.9ECNI11)  , DlStllS i.COPPdSI  ,COR»l  131  ,C3RI  (IS* 

5 PI*I.141*92b535S979 

READ  18, NCC, FREQ, PULE, TEL 
10  FORNATIIS.lc.lS. 01 
RTL* VEL/FREQ 
pull* tel /pule 

10  AFAC*MTL/3bO.O 

QFAC*2.0*PI/RTL 
REF*2.0 
I AOTA* 0 

CALL  CODEIN ( I AOTA 1 
IS  CALL  CORTAL 

DO  400  NNN*1,  NCC 
PRINT  28, FREQ. PULE, TEL, MCC 
20  FORMAT (2bMlCARRIER  FREQUENCY  (H2I  > , 1 P 1 1 1 3 . b , / . 

•29M  »ULSE  FREQUENCY  (dIT/SEC)  = , IP1E 1 5. b, /. 2S 
20  »M  TELOCITY  (NETERS/SEC)  * , 1P1 E1S . b , / . 2 > M NUMBER  OF  SIMULATION1-,  - 

• , 15, / 1 

PRINT  30,(tA(J),J*l.ig>.(IQ(JI,J*l,10) 

30  F ORMAT ( J2M  SHIFT  REGISTER  CONFIGURATION  • , I 01 1 . 2X , 1 u I 1 1 
IBOTA* I AOTA 

25  READ  40.NPN.NCRT, IAOVA 

46  F ORMAT (1I5.5X,1I5,5X,IIS) 

NIOTA* I AOTA -I BOTA 
IFINAUTA.EQ.01  GO  TO  SO 
CALL  SHFTAOTINAOTAI 
30  CALL  CORVAL 

SO  PRINT  bO,  IAOVA,  (iCUI.J-1,  1023) 

bO  FORMATI3SM  RELATIVE  SHIFT  REGISTER  AOVANCF  * ,ll4,//,lCrt  jULJ  CO J- 
•,/.41(lM  ,2515,21) 

PRINT  70, (TCRTL(JI,J*1, 10231 

IS  70  FORMAT (JIHISJLU  COO E CORRELATION  F JNC T I ON , /, OF (1 H ,12113, /II 

RFAO  110,  IPL(JJ),AHP|JJ|,NPI(JJI  .TPIIJJI,  JJM.NPHI 
110  FORMAT  (IE  15.  i , 5)!,  1 FS  . 3 ,S  X , 1 1 5 , S X , 1 7 5 . 3 I 
PRINT  120.NNN.NPM 

120  FORHATT/////,  21M  SIMULATION  NUMjFR  * , 1 IS  , 10  X ,2bHNUM  JE  R OF  -MITTAL 
40  • RATES  * ,lt?,////,7H  NUMBER, 8 X , 1 1 HP AT H LENGTH, 9X» (HAN PL  TTU  I-, 3X, 

•10HNO  OF  INCS,SX,13HI0T  InC  (OEGI./t 
PRINT  131,( JJ ,PL ( JJ) , AMP ( JJJ ,NPI ( JJI , I PI  < JJ> , JJ* 1, NORI 
1 30  FORMAT I IX  , 1 Id  , bX,  1P1E 1 *.  b,  5X,  IP  1E1 5.6,  SX  , 1 1 5 , 1 JX  , OP  IF  9.  3 , 7 I 
A N OR N = 0 • 0 

45  OU  150  11*1, (PR 

IGA((  II)  rPu  (i  1 | / -*ULL 
C H I ( 1 1 1 * P L I 1 1 1 • F l U m I ( 1 j A ( I I | I • J UL  . 

CHlO(II)*CHI( II) 

ANORH* ANuRMAA  MP ( 1 1 1 
S3  158  CONTINUE 

XNORH*REF/(ANORN*204b.OI 
DO  400  11*1 »NPM 
IFIII.EQ.il  CO  T 3 153 
IF (N»I ( I 1 1 .El . 0)  GO  TO  400 
SS  lb 0 r»L!E1*PL(IT) 

CMl  T- N*CMT  ( it  I 
IGTFN*ICAH(II 1 
III«NPI (IIlAl 

IF(npKIII.lQ.u)  GO  TO  170 

bO  OFlANG*TPI(  II I *AFAC/FLOAT(NPI  (H  1 1 

ANSI  = r PI  (I  I)/ FLOAT  (NPI  (III  1 
GO  TO  130 
170  DEL  ANG*  0.0 
ANG 1*0.0 

bS  180  DO  ISO  JJ*1,III 
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70 


76 


*0 


»5 


90 


95 


100 


105 


110 


115 


120 


125 


150 


OJJ. JJ-1 

PHIII=PL,rN»JJ5»OfLANG 
IGANI  1I|=PI  (I  1)  /•’UcL 

Cut (in  a pi i III -float i i cam  i i i i i **011 

OC  ITJ  KK=1,NPK 
CMlO IKKI *CMl ( <<) 

190  CONIIMUF 

ANC  = OJ J*  ANG I 
NPRM=NPN- 1 
OO  200  K*»1 ,NPRM 
K«*<<»1 

OO  200  IL  = K«,N»< 

IFIOMl3<LLI.,t.CMlOIK<ll  Go  TO  2LU 
CHI 1=CHI01LLI 

CHioii.t.1  =c*n  <<<> 
cm  j«ki 

200  CONTINUE 

OC  21b  KK*1,NPH 
AkG*QFAc»Pl  (<KI 
3ESN<hi<|  *STN(  AR3I 

9ECN i A< I * COS  I argi 
210  CONTINUE 

PRINT  220,11. ANG 

220  F0RMATI2  4M1PAKTTAI  a A V E 9<"ING  VAKIEO  = .112,/. 

•374  ANGULAR  PMASt  IwCPEir,!  (JFGRcr.>l  = .IF,.’,//, 

•ROM  CORRELATION  FUNCTION  « IN  RANGt  OF  iNlFpfSIl./l 
PRINT  250 

250  FORMAT  (2«, KLJM  RANGE  AMPLTTUOF  REAL  PA»I  I MAG  PART  ),/| 

NN=1 

NCRVAaO 

OO  290  UK  * 1 .15 
OO  240  LI* 1 , NPM 
NCRXA=NCRVA*1 
COPRlNNI =0.0 
CORI (NNI *0.0 

OELl  = FLOATT<K-l»  ♦ CHIU  ILL  I /HUlL 

OO  250  NM=1,N»H 

RnG***l(  MM  l/PULL  - Ok  LT 

IF  1 RN5.LT.1 .0  IRNG* »NG*  10  25.0 

NUNM=»NG 

CM1  MM=  RNG-FlOA  I (NUMM) 

NUMP*NUMM*1 

IFIMUMP.GT. 10751  NUHP*NUNP-1023 

dET8*AHPINNl*  « I CR  VL  ( NUMM  I ♦ CM  I MM*  C I CR  XL  l NURPI  -ICRXLINJMNI  ll 
GORKI  NNI  *C0RR(NNI  *B£ I B *8  EC NINM I 
COR  I INN»  *C3RI  I NNI  »8k  T b*8E  SN  (MM  I 
250  CONTINUE 

CuRKl NNI  = * (OR M* CORK  I NNI 
COR  1 1 NNI *XN0RH*C0RI«NNI 
OISTI NNI  = PULL • OEL T 

Cu»p(NNI  *SC  RT  ICOMRINnI  *CjRRINNI»CJRI  (NHI»f  OR1  Ir.NI  I 
IFINN.LT. 5)  GO  TO  250 

PRINT  260,  niSKNMI,  COPPINMI  .CORRCIHI  ,COPI  (NHI  ,N«=  1.  T) 

260  FORMATim  , T(  1P%:  10,2, 5X1  I 
NN=  0 

200  NN*NN ♦ 1 

IF(N'.RNA-*lC4V  .GE.  5IGU  TO  795 
290  CONTINUE 

TFINCRXA.GF.Ni.RlM  GJ  TO  59C 
295  PRINT  500 

500  FORMAT!///, 5IM  CORRELATION  FUnCTIJN  IUUTGIl""  F ~ I C~  I ^L « ^ ' N 

•IERESII ,/l 
PRINT  250 
Do  510  NN* 1,15 
roRP<NNI*0.9 
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13S 


140 


14P 


IS  u 


IS, 


15  j 


1 n -> 


i /j 


1 25 


iso 


180 


COMRINNKd.O 
CORK  NNl *0.3 
310  CONTI N UE 

UO  JjO  KK=1o. 1023 
00  3 > 0 cL*l.NPM 
NCKOA; N jPVA  »1 
CPRKA=0. 0 
COMI4*0.1 

0ELI  = FL0AH«-11»CHIj(lL12PJLI. 

00  J?0  H’UliSi'il 
9HS* ’Ll MM1 /BULL -OtL t 
nUHMKNG 

CMINN*9NG-rLJAT INUNMi 
IKnjMn.LT.il  NjMM  = nuim«i:j3 
NU1P*NUM"M 

IFINJMP. GT. 10231 NO HP  = N UN P-1 1 23 

oTTil:  4 HP  (Mil*  I ICMXLl noHMI  ♦ CnlPMKIC’VLI  l.UM'l  - 1%  jp»  I » 1 

00><-?4*C0Po4»lLr3*3eCNCMMI 
:OkI**„o«>U*3U  l3»atSN(.1“l 
32n  CONTINUE 

«,J<MA  = XNjF"»C')-m 
C0°i4:XNjOM»J.  0014 
PI  G I 4 = POLL  * Or  L T 

CORPA  * ,»QRT  ( CO  9PA  • CORK**  *C  J1,  I A *C0N1 4 I 
)j  3 4V  MM-  l , l a 

if«Cj<°a.u .; jspimmi i&o  fo  is, 

IF  IIM.&C.I  .1  r'Z  TO  1-. 

NM* 1 i -MM 

U.  3 Jj  I.N*1  , (“ 

N4-1  j-.,N 
in-  1 S-NN 

CoMP I3II>C JP  ' ( M 4 | 

OLMMt  1 II  *r.?-U  [NAI 
COM  1 1(4  il  .-CJM  ( NAI 

0 1 G T (Nril  rnTOT  ((4  4) 

33P  cjnrtsjr 

34"'  CPMPIM“>sr^3<3 
CC9HMP1  Jf:’.  4 
COM'-1  (MmI  rfOMJ  4 

nioi  (mmi  roi  ;i  a 

IK  4-<V4-N'»i/.«l.01  r.u  TO  3oO 

so  n iso 

3 4 S CuuTlNUE 
.35  0 CONTINUE 
360  NNAA*LLMK<-tbl*NPM 

IF 1NNAA.GF.15INNAA*15 

PPiNf  3 7 3,  niSKNNI.COPPlNNI  .C09RINN1  ,C09I  (NNI  ,N  ,=  l.NNAAl 
3 7 C FOPMAI (3(1<.1D4E10. 2,2X1! 

380  CONTINUE 

PsIIIKPLTFN 
CHl(III*CHIT£N 
IGAM(  IIMTOTFN 
400  CONTINUE 
STOP 
ENO 


1 SUBROUTINE  COPi/Al 

COMMON  IC(1  02  II  • ICMLC 10231  , IA(10I  , IB  1 101 
DO  100  KOIF*1,102J 
CALL  COOEC39I KDIF.ICRVLCKOIFll 
3 100  CONTINUE 

RETU»N 
END 


■NT 


34 


I 


1 


5 


10 


19 


20 


25 


I 


10 


35 


50 


%s 


CUJNJIJIINF  33  OF  I N C NA  0 W A I 
DIHFNSI3N  I JNM  111  , I5R  I ( 1 01 

cohn-in  iccio23i,icrwiciq23i,iacigi,ibcioi 

»fll)  lQO.CTACJI  « J *1 , 1 0 1 « 1 1 J ( J I , J = 1 , 1 0 1 
10"  FU*-C*iA  T 11011*1  OX.lWllI 

on  no  j=i,io 

1 SR  A ( Jl  = 1 
ig»bi ji « i 

150  CONTINUE 

ENT3T  SMFIA3V 

200  IFINAOWA.&r.O)  S3  to  211 
TF CNAJWA.ti.OI03  TO  250 
NAQWA=1Q23*NA3WA 
GO  f 3 200 

210  00  200  J=l,NAOWA 
NA*0 

00  220  <=1,10 

NA  = NA»1A  C<l M SRA IK| 

220  CONTINUE 

00  2J0  <=1.3 
C * 1 1 - < 

LL*10-« 

Ij<ACL»=I59ACLlI 
230  CONTINUE 

1 SR  A 1 1 1 «i«A-2*  I NA/2  I 
200  CONI  I NOE 

250  00  500  J=  1 , 10  2 J 
NA=  0 
N8*0 

00  300  <=1,10 
NA=NA. I A C <1  « I SNA| <) 

NB  = N3*TBKI»ISR3C<1 
300  CONTINUE 
ICCJ)*-1 

IFCISRAC10I.lC.I5R3(1QIIICCJI=1 

00  A 0 0 <>1,3 
L*ll-< 

LL= 1 !•< 

1 SR  A C L l = IS®  A C ELI 
ISR3 ( L I = IS»3CLL» 

400  CONTINUE 

ISRA < 11 »NA-2»  CNA/21 
ISRB(ll=NB-2»  CNB/21 
500  CONTINUE 
RETURN 
END 
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5 


10 


15 


SUJR3UT  1NL  GJULCJNKOAF,  IS  UNI 

CONNON  ICC102  51  , luRWLC  10211  ,IA  (101  , I B C 1 0 I 

<OIF  = <OAF 

10  IFKHF.wF.ul  Co  TO  20 
<OIF=<OIF*10?3 

go  n io 

20  IFC<JIF.&1.1023|<OIF=<jIF-1J2J»(kuIF/1023) 
I 5UN  = 0 

00  100  1=1,1323 
< = I ♦<  OIF 

IF(<.GT.102IX»<-102i 
I SUH« l SUN* T Cl  1 1 • ir (<l 
100  CONTINUE 
RETURN 
FNO 
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